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ABSTRACT 

Disk galaxies at high redshift have been predicted to maintain high gas surface densities due to 
continuous feeding by intense cold streams leading to violent gravitational instability, transient features 
and giant clumps. Gravitational torques between the perturbations drive angular momentum out and 
mass in, and the inflow provides the energy for keeping strong turbulence. We use analytic estimates 
of the inflow for a self- regulated unstable disk at a Toomre stability parameter Q ~ 1, and isolated 
galaxy simulations capable of resolving the nuclear inflow down to the central parsec. We predict an 
average inflow rate ^ 10 M© yr~^ through the disk of a 10^^ Mq galaxy, with conditions representative 
of z ~ 2 stream- fed disks. The inflow rate scales with disk mass and (1 -I- z)^^^. It includes clump 
migration and inflow of the smoother component, valid even if clumps disrupt. This inflow grows 
the bulge, while only a fraction ^10"'^ of it needs to accrete onto a central black hole (BH), in 
order to obey the observed BH-bulge relation. A galaxy of lO^^M© at 2 ^ 2 is expected to host a 
BH of ^IO^Mq, accreting on average with moderate sub-Eddington luminosity ix'^lO^^"^'^ ergs~^, 
accompanied by brighter episodes when dense clumps coalesce. We note that in rare massive galaxies 
at z r^6, the same process may feed ~ 10^ M© BH at the Eddington rate. High central gas column 
densities can severely obscure AGN in high-redshift disks, possibly hindering their detection in deep 
X-ray surveys. 

Subject headings: galaxies: formation — galaxies: nuclei — galaxies: active — galaxies: high-redshift 
— X-rays: galaxies 



1. INTRODUCTION 

The growth of massive black holes (BH) and the as- 
sociated Active Galactic Nuclei (AGN) are commonly 
assumed to be driven by gas inflows from gal axy mergers 
(jDi Matteo et all 120051: iHopkins et aIll2006D . However, 
baryons are fed into high-redshift galaxies through cold 
streams, in which the contribution of smooth gas and 
small clumps is larger than that of mergers of mass ra- 
tio greater than 1:10 (jDekel et al.l I2009at iBrooks et all 



[2OOI . 

Cold accretion and high gas fractions in high-redshift 
disks trigger a violent instability, often characterized 
by giant star-forming clumps. Using analytic estimates 
(Section 2) and simulations (Section 3), we address the 
possibility that the gas supply for BH growth is provided 
by instability-driven inflows, and evaluate the obscura- 
tion of AGN in this model (Section 4). 



2. INFLOW IN HIGH-REDSHIFT UNSTABLE DISKS 
Many z ^1 galaxies of stellar masses w 10^"^^"'^ Mp 







are rotating disks with g iant clumps of ^ly ^M0 each 
(jElmegreen et alj 120071: iGenzel et all I2011D . In most 
cases, o n-going mergers a re not favored by gas kine- 
matics (IGenzel et all I2008D and photometric properties 



(lElmegreen et alJ 2009D. Instead, high gas fr actions of 
50% (jDaddi et alJ l2010l : iTacconi et al.ll2010l) and high 



gas velocity dispersions (jForster Schreiber et alJ l2009( ) 
are consistent with gravitational instability and self reg- 
ulation at Q ~ 1 . 

The occurrence and persistence of gravitational disk 
instability at high redshift is a natural result of the 
high surface density of the cold disk component, which 
is maintained by the continuous gas su pply through 
cold cosmic streams (jDekel et alJ l2009a| ) . The insta- 
bility is predict ed to be self regula ted in steady state 
for a few Gyr (jDekel et alJ I2009H DSC09), and cos- 
mological simulations reveal clumpy disks that resem- 
ble th e observed ones (|Agertz et al.ll2009t iCeverino et alJ 
I2010D . If giant clumps survive stellar feedback for a 
few 1 0^ yr, th ey migrate in to a central bulge (|Noguchil 
119991 : iBournau d et al.ll2007D . as supported by observa- 
tions (jGcnzcl_£^_aI. J2008 ^^_Wpica]_chm to sur- 
vive (jKrumholz fc DekellboialGenzel et al.ll2011h . while 
the most luminous clumps could undergo faster disrup- 
tion ([jylurray ct al. 2010: Gcnzcl ct al. 20Tll). Neverthe- 
less, the inflow is a robust feature of the disk instability, 
driven by gravity torque s independent ly of the surviv- 
ability of bound clumps (jGammiell200ll DSC09). 

A key element in the self- regulated instability at Q^l, 
given the high surface density, is maintaining strong tur- 
bulence, with a velocity dispersion ~50km s~^for circu- 
lar velocity ~ 200km s~^, as observed in high-z disks. 
Supersonic tur bulence dissipat es in a few disk vertical 
crossing times (jMac Lowlll999f ). and should therefore be 
continuously powered. An obvious source is the gravita- 
tional energ y released when mass flows toward the cen- 
ter (DSC09;MmeireenX£urE3[20i3). Assuming that 
the incoming streams deposit mass and energy at the out- 
skirts of the disk where they do not contribute directly 
to driving disk turbulence, as is suggested by cosmologl - 
cal simulations (|Agertz et al.ll2009l:ICeverino et al.ll2010[ ). 
the main energy source is the inflow within the disk 
(|Krumholz fc BurkertI [20101 ). This energy is transferred 
to velocity dispersion by gravitational torques involving 
massive clumps, transient features and smoother regions 
of the disk. Stellar feedback is unlikely to be the main 
driver of the turbulence, e.g., because the velocity disper- 
sion is observed not to be t ightly correlated with proxim- 
ity to star- forming clumps (jForster Schreiber et al.l[2009l : 
IGenzel et aT1[20Tlh . 



The baryonic inflow rate Mb is estimated to be: 



td W 



(1) 



where Md is the disk mass, td is the crossing time, V 
is the circular velocity and a is the ID velocity disper- 
sion. At this rate, the power released by the inflow down 
the potential gradient, M^V'^, compensates for the dis- 
sipativc loss rate, Afgcr^/tdis, where Mg ~ O.SAfd is the 



disk gas mass and idis ~ 3id (jCacciato et al.|[201lD El- 
Similar estimates are obtained from the rate of energy 
exchange by clump encounters (DSC09, eqs. 21 and 7), 
and from the angular-momentum excha nge among the 
transient pe rturbations in a viscous disk (jGammiell2001l : 
IGenzel etal.- 2008. DSC09, eq. 24). Adopting a/V^O.2 
and td~30 Myr at z = 2, the inflow rate is 



Ah ^ 25Afb,ii(l 



\3/2 



(2) 



where A/b.ii is the disk mass in 10 Mq, and (1 + z)^ — 

(1 + z)/3. Note that the (1 + z)^ scaling corresponds 
to the disk crossing time, proportional to the halo cross- 
ing time at fixed overdensity in a matter-dominated Uni- 
verse. 

A continuity equation for the disk gas, being drained 
into star formation, outflows and bulge and replenished 
by streams, yields at z < 3 a steady state with A/b about 
one third of the cosniological baryonic accretion rate 
(|Bouche et a l. 2010; Krumholz & Dekel 20^1) ■ With the 

latter being - 80 M12 (1 + z)I^^Mq yr"! (jNeistein et al.l 
120061 Afi2 being the halo mass in IO^^Mq), this esti- 
mate is consistent with Eq. ^ at z = 3, with a different 
redshift dependence for the cosmological inflow rate. 

3. SIMULATIONS 

We perform six idealized simulations using t he Adap- 
tive M esh Refinement (AMR) code RAMSES (jTevssied 
l2002f l. with a density-dependent grid refinement strategy 
and a barotropic cooli ng model producing re alistic phase- 
space ISM structure (jBournaud et al.ir2010[ ). The initial 
parameters of the six simulations are listed in Table |TJ 
The smallest cell size is eAMR=1.7 or 3.4 pc in high- 
resolution (HR) runs and 10.2 pc in medium-resolution 
(MR) runs. Cells are refined if they contain gas mass 
larger than 2 x 10^ Mq, or more than 20 particles, or if 
the Jeans length is resolved by less than four cells. The 
initial number of dark-matter and stellar particles are 
4 X 10^ and 6 x 10^ in HR and MR, respectively. Star 
formation occurs above a density threshold of 5 x 10^ 
and 4 x 10'^ cm~'^ for HR and MR, with a star-formation 
efficiency per local free-fall time of 4%. A fraction of the 
supernovae energy rysN = 0.5 is released as kinetic feed- 
back. For MR4, feedback is boosted up after 250 Myr to 
VSN = 5 in order to disrupt the clumps. 

The simulations start as pre-formed disks in dark mat- 
ter halos. The exponential disk is truncated at i?d, 
with scale-length i?d/4 and i?d/l-5 for stars and gas, re- 
spectively. The dark-matter halo has a Burkcrt profile 
with a scale-length i?d/3. The dark-matter mass frac- 
tion within i?d is 0.25. The key parameter distinguish- 
ing high-redshift disks is the gas fraction, set to 50%, as 
observed at z ~ 2. A small bulge containing a fraction 
/bulge of the stellar mass is assumed. 

Within one orbital time, the violent gravitational insta- 
bility, with Q <1 in the initial conditions, develops giant 
baryonic clumps and supersonic turbulence. The simu- 
lations resolve dense gas structures up to ^ 10^ cm^"^, 

^ we used 0.5/3 Si 0.2 in Eq. 1 



TABLE 1 

Simulation parameters, average gas inflow rates Mg between t = 250 and 400 Myr, and predicted value /gA/b from 

Equation [21 in Meyr^^. 
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Fig. 1. — Surface density of gas and stars at two snapshots of simulation HRl. The gas is shown face-on (top) and edge-on (middle), the 
stars are shown face-on (bottom). To mimic optical images, the luminosity of each stellar particle is constant during its first 10 Myr then 
decreases as t~^'^. Clumps Ai and A2 merge together between the snapshots. 
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Fig. 2. — Relative change in angular momentum per rotation 
period at radius r in the disk of run HRl, as evaluated from the 
gravitational torques. 

with a turbulence cascade in the ISM probed by the spec- 
tral analysis of model HRl (jBournaud et al.ll2010L Sec- 
tion 3.4 and Fig. 15). These isolated models allow us 
to follow the evolution of unstable disks for a couple of 
orbital times, i.e. a few 10^ yr, after which external gas 
supply should be included. 

These simulations demonstrate that the instability- 
driven inflow is a generic feature, not limited to the mi- 
gration of clumps. Figure [1] shows two snapshots of run 
HRl. At t = 210 Myr, giant clumps have formed in the 
outer disk and the turbulent speed has reached a steady 
state with a mass-weighted average cr ~ 65 km s~^. 
About half of the gas mass is located between the clumps, 
the central kpc containing initially only 3x10^ Mq of gas. 
At t = 330 Myr, the giant clumps have barely started 
their inward migration but the central kpc already con- 
tains 1.5 X 1O^M0 of originally diffuse gas that was driven 
in by torques. 

Figure [2] shows the relative change in angular mo- 
mentum per rotation period due to the gravity torques, 
AL/L = 27rFt/Fr, where Ft and F^ are the tangential 
and radial components of the force, averaged in radial 
bins. In the inner 4 kpc, most angular momentum is re- 
moved within one rotation, demonstrating the efficiency 
torquing driving the inflow. 

Table 1 lists the gas inflow rates through spherical 
shells at r = l kpc and 150 pc, averaged between t = 250 
and 450 Myr, i.e. between the onset of clump forma- 
tion and their central coalescence. The simulated gas 
inflow rate at 1 kpc is in the ballpark of the theoreti- 
cal estimate of Eq. [5J smaller by up to a factor of two, 
probably because the replenishment by fresh gas is ig- 
nored in the simulations. The inflow rate is not af- 
fected by the stronger feedback and clum p disruption 
in run MR 4. Indeed, in the simulation by IGenel et al.l 
(J2010D in which the clumps are disrupted by enhanced 
outflows, the inflow rate is similar to that in their ref- 
erence simulation where the clumps survive, reaching 
in both cases a peak gas surface density higher than 
10^ Mq pc-2 pa 1023cm-2. 

Figure |3] shows the net gas inflow rate across spheri- 
cal shells about the stellar center of mass. Local inflow 
rates exceeding 10 M0 yr~^ are seen beyond the inner 
1.5 kpc, dominated by clump migration. The inner kpc 
is dominated by smooth inflow at a rate of ~ 5Mq yr~^. 
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Fig. 3. — Radial profiles of the gas inflow rate. Top: not 
inflow (inflow minus outflow) through spherical boundaries for run 
HRl at three different times and for run MRl. It is dominated 
by clump migration at large radii (note that clump B is moving 
outward), and by smooth inflow of 1 — TMQyr"-*^ between 70 pc 
and 2 kpc. Below ~ 30 — 50 pc the inflow rate drops, but it remains 
at 1 — 2M0yr~^ in the inner 30 pc. Bottom: Net inflow rate 
through cylindrical and spherical boundaries for HRl, comparing 
inflow and outflow. The drop in the net spherical inflow rate at 
small radii is dominated by the geometrical effect, and outflows 
further reduce the net inflow. 

The mild decline between 1 kpc and 50 pc can be partly 
explained by a geometry effect, as the sphere radius be- 
comes smaller than the disk thickness. This is demon- 
strated by the inflow rate through cylindrical boundaries 
being rather constant with radius (Fig. 3, bottom panel). 
A large fraction of the inflow takes place above and be- 
low the mid-plane in the thick gaseous disk. An inflow 
rate of 1 — 2M0yr~^ nevertheless persists through the 
central few parsecs across spherical boundaries. We no- 
tice that non-negligible outflows, presumably from stellar 
feedback, further reduce the net inflow rate in the central 
tens of pc. A comparison of runs HRl and MRl indicate 
no resolution effects on scales of 10-100 pc. 

This gas inflow fuels star formation in the bulge. 
Galaxy HRl at t = 300 Myr has a total SFR of 
132 Mq yr^^, of which 16 M© yr~^ is in the central 
kpc. This is less than the SFR of 42 and 31 Mq yr'^ 
in clumps A and B, but comparable to the 14 Mq yr~^ 
in clump C. At t—450 Myr, the stellar bulge mass has 
grown from 9.8 x 10^ Mq to 2.2 x IO^^Mq. The persis- 
tence of star formation in the bulge as a result of disk 
instability is consistent with observations showing that 
bulges in clumpy disks are bluer than bulges in smooth 
disks (|Elmegreen et aLl f" 



4. AGN OBSCURATION IN GAS-RICH DISKS 

The central BH resides in a dense and thick gaseous 
disk (Fig. [1]), which could obscure an AGN. To quan- 
tify this, we computed the column densities over 60 ran- 
dom lines of sight through the centers of galaxy HRl 
at t=330 Myr and HR2 at t=380 Myr. The results are 
insensitive to the choice of center, selected to be either 
the center of mass of old stars or the point of maximum 
vorticity. The foreground hydrogen column density iVjj 
was estimated over a cross-section corresponding to the 
minimum Jeans length, 4 x camr- 

Figure |4] shows the distribution of A^h for runs HRl 
and HR2 and for various sets of inclinations. 15% of the 
lines of sight reach Compton thickness, log(iVH) > 24.1, 
and 50% have log(A^H) > 23, i.e., strong obscuration. 
The obscuration tends to be higher in the more massive 
disk HRl , where Compton thickness can be reached even 
for face-on orientations owing to the ~ 1 kpc thickness of 
the gas disk. The lower- mass galaxy HR2, perhaps more 
representative at z < 2, can reach Compton thickness 
almost exclusively for edge-on orientations. 

If most gas in the central parsec lies in a torus, un- 
resolved in our models, the obscuration on edge-on pro- 
jections may be further enhanced and the dependence 
of obscuration on inclination may be stronger. Signifi- 
cant obscuration would still be present along most lines 
of sight, because the column densities measured in the 
simulations are on average not dominated by the central 
10 pc. 

5. DISCUSSION: ACCRETION ONTO THE BLACK HOLE 

The isolated simulations presented here, using high 
gas fraction representative of z -^ 2 disks, have a disk 
inflow rate consistent with cosm ological simulations o f 
high-redshift stream- fed galaxies (jCeverino et al .11201 Of ). 
They reach resolutions better than 2 pc, hence resolving 
the nuclear inflow down to scales at which other processe s 
drive the actual small-scale BH accretion ()Combesll2001[ ). 
Only a fraction of the inflowing gas needs to make it all 
the way into the BH. The issues are similar to the case 
of merger-driven fueling, namely getting rid of angular 
momentum while avoiding excessive star formation and 
AGN feedback. This is beyond the scope of our paper, 
and we limit ourselves to heuristic estimates. 

One can assume that the local relation between the 
BH mass Mbh and the bulge properties (mass Mbig and 
velocity dispersion a) is crudely valid at high redshift. If 
we adopt MBn/A^big ~ 10"'^ at z = 0, assume it scales 
as a^ and allow a cosmological scaling cr^ oc (1 H- z), 
then we obtain MBn/Mbig -- 3 x 10"^ at z -- 2. A 
z ^ 2 galaxy of baryonic mass 10^"'^ M© hosts a BH of 
^ 10* Mq, while in a gravitationally unstable steady state 
it is typically half bulge and half disk (DSC09). Ac- 
cording to eq. ©, the average inflow rate through the 

disk into the bulge is Mbig ~ 12.5M0yr~^. If the ra- 
tio of average accretion rates into bulge and BH follow 
the ratio of the corresponding masses, the BH accretes 
on average at A/bh ^ O.O4M0 yr"-'^, which is 1-2% of 
the Eddington rate (assuming a 0.1 efficiency for mass- 
to-energy conversion). The corresponding bolometric lu- 
minosity is 2 X lO'^^ergs"^. With typically 1-5% in X- 
rays, we estimate on average Lx ~ 10^^~'^'^ergs~^, scal- 



ing with galaxy mass and with (1 -|- z)^'^. While the 
average luminosity would be modest, short episodes of 
higher accretion rate, possibly up to the Eddington level, 
occur during the central coalescence of migrating giant 
clumps - which could a lso bring with them seed BHs 
(jElmegreen et al.ll2008ah . 

At very high redshift, the same pr ocess could feed 
brigh ter AGN in rare massive systems (|Di Matteo et al.l 
120 lit) , provided that unstable disks indeed form there. 
The average inflow rate increases with redshift (eq. 2), 
such that a z = 6 — 10 disk can support continuous accre- 
tion at the Eddington rate in self-regulated quasar mode 
if A/sH/A/d is --0.04 - 0.08, for the same - 10~^ fraction 
of the inflow rate accreting onto the BH, and assuming 
outflows would not affect the cold gas streaming inward. 
Then a seed BH of 10^ M© at z ~ 10 can grow expo- 
nentially by 10 e- folds to 2 x 10^ Mq at z~6, possibly 
explaini ng very massive BHs growing in gas-rich systems 
at z > 5 (jMortlocket al.li2011 : Treister et al. 2011) . The 
typical baryonic mass of the host galaxy is ~ 10^^ M©, 
and if the halo mass is a few times 10^^ Mq, the comov- 
ing number density of such bright z > 5 AGNs would be 
-10-** to 10-6 Mpe-^ 

6. CONCLUSION 

Violent gravitational instability in high-redshift disk 
galaxies naturally leads to the growth of a bulge and 
a central BH; it can both feed an AGN and obscure it. 
This stems from the developing picture where many star- 
forming galaxies at high redshift are rotating disks with 
high gas fractions continuously fed by cold streams. Such 
disks develop instability that involves transient pertur- 
bations and giant clumps. Gravitational torquing among 
these perturbations lead to mass inflow, which provides 
energy for maintaining the strong turbulence required for 
self-regulated instability at Q ~ 1 . The cosmological in- 
flow rate along cold streams sets an upper limit, but the 
actual inflow rate toward the nucleus is determined by 
the disk instability on galactic scales. 

High-redshift disks are very different from nearby spi- 
rals, the key difference being the much higher gas frac- 
tion. The instability then operates on short dynami- 
cal timescales, and drives an intense inflow through the 

disk at the level of - IOMq yr-iMb,ii(l + zfj'^ . In 
low-redshift disks, the instability is a secular process 
limited to non-axisymmetric modes with considerably 
weaker torquing: bars are invoked in some AG N models 
(jFanidakis et al.l 120111: iBegelman et al.l I2006D but need 
about ten rotation pe riods to convey some gas inward 
(jBournaud et al.ll2005[ ). Gravitational torquing in high- 
z disks is much more efficient and involves richer gas 
reservoirs. 

Disk instability at z — 2 can thus funnel half of the 
disk gas toward the center in 2 Gyr. This is similar 
to th e mass inflow in a major merger (jHopkins et al.l 
120061 ). but spread over a ten times longer period, result- 
ing in a lower average AGN luminosity, with higher duty 
cycle, and high obscuration. This process could domi- 
nate, as the gas mass in the relatively smooth cosmolog- 
ical streams is larger than that associated with mergers 
- less th an 10% of the cosmological infall is in major 
mergers (jDekel et al.„2009a; .L'Huillier et al.„201ll ). The 
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Fig. 4. — Distribution of gas column density along random lines of sight in front of a central AGN. We show separately the distributions 
for face-on (i < 30°), intermediate (30° < i < 60°) and edge-on (i > 60°) orientations, the total fraction in each of these sub-samples being 
normalized to 33% for clarity. 

main prediction is thus that many high-z AGNs should 
be hosted by star-forming disk galaxies, composed of 
clumpy disks and growing spheroids. Many would be 
moderately luminous, detectable by Chandra deep sur- 
veys a.t z — 2 ii un-obscured, but the high expected ob- 
scuration makes their detectability dubious, especially 
for galaxies belo w 10^^M(7). These BHs can grow inside 
classical bulges ([Elmegreen et al.ll2008b[ ). 

Observationally, the fraction of X-ray-detec table AGN 
indee d declines with decreasing stellar mass (jXue et alj 
I2010D . Many AGN are hig hly absorbed, indi vidually un- 
detectable in X-rays (e.g.. lDaddi et al.ll2007D . and many 
of these obscured AGN have modera te intrinsic luminos i- 
ties and lie in star- forming galaxies (iJuneau et al.|[20Tll ). 
Host galaxies of moderate-luminosity AGNs do not show 
an excess of major- merger signatures (jGrogin et al.ll2005t 
IGabor et al.l 120091) and often have disky morphologies 
([Schawinski et alTl2011|) . Small mergers may also help 
feedin g the BH while leaving the disk intact ([Shankarl 
|2010[) . but they are a component of the cosmological 
streams, and can be considered to be part of the inflow 
within the disk. 

The violent instability phase should end when the cos- 
mological accretion rate declines and the system becomes 



stellar dominated. We expect less massive galaxies to 
remain unstable for longer times, because they retain 
higher gas fractions — one o f the aspects of the down- 
sizing of star formation (e.g.. IJuneau et ai]|2005[) . This 
can result from the regulation of gas consumption in 
smaller galaxies (|Dekel fc Silklll986t iKrumholz fc Dekell 
120 lit) and from the continuation of cold accretion to 
lower redshift for lower-mass halos (jDekel fc BirnboimI 
|2006|) . It induces an inver se gradient of gas fraction with 
galaxy mass (as observed. IKannappanll2004l ) and a down- 
sizing in gravitational instability which could result in 
downsizing of BH growth. This is a longer growth phase 
into later redshifts in lower-mass galaxies, with ^ IO^^Mq 
galaxies growing BHs mostly at z > 1 in their unstable 
disk phase, while clumpy disks of ~ 1O^°M0 may show 
AGN activity even after z^l. 
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